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Abstract 

The macroscopic behavior of dissipative stochastic partial differ- 
ential equations usually can be described by a finite dimensional sys- 
tem. This article proves that a macroscopic reduced model may be 
constructed for stochastic reaction-diffusion equations with cubic non- 
linearity by artificial separating the system into two distinct slow-fast 
time parts. An averaging method and a deviation estimate show that 
the macroscopic reduced model should be a stochastic ordinary equa- 
tion which includes the random effect transmitted from the micro- 
scopic timescale due to the nonlinear interaction. Numerical simula- 
tions of an example stochastic heat equation confirms the predictions 
of this stochastic modelling theory. This theory empowers us to better 
model the long time dynamics of complex stochastic systems. 

Key Words: Stochastic react ion- diffusion equations averaging, tight- 
ness, martingale. 
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1 Introduction 

Stochastic partial differential equations (spdes) are widely studied in model- 
ing, analyzing, simulating and predicting complex phenomena in many fields 
of nonlinear science [8], [9l [191 e.g.]. 
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Reaction-diffusion equations (rdes) are important mathematic models 
naturally applied in chemistry, biology, geology, physics and ecology etc. 
Such equations can be obtained from microscopic particle systems under the 
so called hydrodynamic space-time scaling limit [5] . When taking fluctua- 
tions under the hydrodynamic limit into account, an additive noise appears as 
correction term to the reaction-diffusion equation [I]. This reaction-diffusion 
equation with noise is viewed as an equation describing intermediate level 
between the macroscopic and microscopic ones. A deterministic equation 
appears in the macroscale only if the randomness averages out. 

However, for nonlinear complex systems, random effects may survive av- 
eraging and thus be fed into the macroscopic system [T^. Especially when 
we consider the macroscopic behavior of solution on a long time scale, such 
random effects should not be neglected [121 [151 El e.g.]. Furthermore, macro- 
scopic turbulence may need to be modeled as noise in, for example, geophys- 
ical fluid dynamics. 

Blomker et al. [21 HI e.g.] recently studied amplitude equations for spdes 
with cubic nonlinearity, which is proved to be a stochastic Landau equation. 
But in the amplitude equation the fluctuation from the fast modes disap- 
pears when the noise just acts on fast modes. For spdes with quadratic 
nonlinearity, Roberts [13] derived the amplitude equation for one dimen- 
sional stochastic Burgers equations by calculating a stochastic normal form 
model. Blomker et al. [3] also gave a rigorous proof for a general spdes 
with quadratic nonlinearity by a multiscale analysis; they showed that the 
amplitude equations include the fluctuation from the fast mode due to the 
nonlinearity interaction. 

This paper considers react ion- diffusion systems driven by a noise which 
is homogeneous in space and white in time. Here we are concerned with 
the dynamics of the system on a long time scale. For this, a scale transfor- 
mation separates the system into slow and fast modes. Then we derive a 
low dimension macroscopic system which provides the long term dynamics. 
And the low dimensional macroscopic system includes a noise term which is 
transmitted from the fast modes due to nonlinear interaction. 

For definiteness, let the non-dimensional / = (0, vr) and L'^{I) be the 
Lebesgue space of square integrable real valued functions on / . Consider the 
following non-dimensional react ion- diffusion equation 

dtW = dxxW + f{w) + adtW on < x < vr , (1) 
w = on X = , vr , (2) 

where f{w) represents a nonlinear reaction and W is an L'^{I) valued Q- 
Wiener process defined on complete probability space [Q, J-', P) which is de- 
tailed in the next section. Our aim is to study the behavior of solutions 
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to ([I])-® over a large timescale, say e^^ for small e > . For some fixed 
integer > , we split the field w into N low wavenumber modes and the 
remaining high wavenumber modes. For this non-dinesional problem, denote 
by {ek{x)}k = {sm{kx)}k the orthonormal eigenvectors of d^x- Define the 
projection operators onto 'slow' and 'fast' modes respectively 

N 

Pat- = ^efc(x)(efe(x) , ■) , Qj^=X-Vn, (3) 
fc=i 

where I is the identity operator on . Then defining u = Vnw and 

V = Qnw , the RDE (P) is identical to the following coupled equations 

dtu = d,,u + VNfiu + v) + aVNdtW, (4) 
dtv = d,,v+QNf{u + v) + aQNdtW, (5) 

whence w = u + v . In order to completely separate the time scales we modify 
the above system by introducing a 'high-pass filter' Aj^ defined by 

An = - (1 - t)VNdxx = {Qn + (^VN)d^x ■ 

Observed that when e = 1 , the physical case, the operator = d^x as 
appears in ([T]) and dl])-®; but when e = 0, the operator Ajq is a pure 
high-pass filter with null space spanned of the 'slow' modes. 

For the moment also assume the noise acts only on the high wavenumber 
modes, that is VjyW = . And modify the system dl])-© to 

dtu' = ANu' + PNfiu' + v'), (6) 
dtv' = ANv' + QNf{u' + v') + aV^QNdtW . (7) 

Here the choice of ^/e, the factor in front of noise term, ensures the fast 
modes, solutions of (JTj), remain of order 1 as t ^ oo and e — > for any fixed 
. Note that when e = 1 , ©-(CD is identical to O)-© and 1^. We aim to 
use analysis based upon small e to access a useful approximation at e = 1 . 

Section [5.11 proves that, for small enough e > 0, high modes v'^{t) is 
approximated by ^/erj^{t) over long timescales (1/e) where ?7* is the stationary 
solution solving the linear stochastic partial differential equation 

dtT] = AnT] + aQNdtW . (8) 

Consequently our careful averaging proves that the macroscopic behavior 
of u''{t) is described to a first approximation by y/euN{^t) which solves the 
following finite dimensional, deterministic system 



dtUN = dxxUN + VNkiuN) 



(9) 
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Here, the average 

- 1 /"* 

/o(MAf) = lim - / fo{uN + r]*{s))ds , (10) 

and /o is the cubic component in /, see detail in next section. Usually © 
is called the averaged equation and in probability the following convergence 
holds (Section l5.2p : for any T > there is positive constant C 

sup \u^{t) - y^UN{et)\L2(D) < Ce (11) 
0<t<e-iT 

for small e under some proper conditions on the initial value. Further, the 
martingale approach of Section 15.31 shows that small Gaussian fluctuations 
generally appear in these slow modes on the timescale e~^. The approach 
shows that fluctuations are transmitted from the fast modes by nonlinear 
interactions. The last Section [6] confirms these theoretical predictions by 
comparing them to numerical simulations of a specific stochastic reaction 
diffusion equation. 



2 Preliminaries and main results 

This section gives some preliminaries and states the main result. First we 
give some functional background and some assumptions. 

Let H = LF'{I)- Denote by A the second order operator dxx with Dirichlet 
boundary on / and let {ei}°^i be a eigen-basis of H such that 

Aci = -atCi , z = 1,2, . . . , 

with < «! < a;2 < 0^3 < ■ ■ ■ . For any 5 > , introduce the space Hq = 
D^A^^"^), which is compactly embedding into H. And let H^^ denote the 
dual space of Hq. The usual norm defined on is written as || • H^. And 
for 6 = and 1, the corresponding norms are written as | ■ | and || ■ || 
respectively. Denote by (■, ■) an inner product in H such as the inner product 
{u, v) = (2/7r) uv dx . And for positive integer A^, denote by Hn the space 
spanned by {ei, . . . , cat} and by H]i^ the space spanned by {e^+i, . . .} . 

We are given a complete probability space (f2, JF, {J-'t}t>o, P)- Assume W 
is an if-valued Q-Wiener process with operator Q that commutes with Vn 
and satisfies 

A^c^ , z 1,2,..., 
with Ai = , 2 = 1, 2, . . . , A^ and Ai > , 2 = A^ + 1, . . . . Then 

oo 

w{t)= V^A{ty^ (12) 

i=N+l 
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where component noises 



P,(t)^^{Wit),ei), i^N + 1,..., 



are real valued Brownian motions mutually independent on (Q, P). More- 
over we assume 

oo 

trg = ^ Aj < oo. (13) 

i=N+l 

In the following we denote by Q2 = {I ^ 'Pn)Q = Q ■ For the nonlinear 
term / we assume the following hypotheses (H) 

1. / : M — > R is smooth and has the following form 

/(0 = /o(0+/i(0, eeM, 

where foiO = ~Co^^ , Co is some positive constant, /{"(O) = and 

/i(oe<o,eeM. 

2- /'(O < for some positive constant Ci . 

3. There are positive constants C2 , C3 and positive integer p such that for 

1/(01 <c2ier^-^+c3, fm<c2e'+cs. 

Then define 

m = ^/(^o (14) 

which is well-defined as e — > for any ^ e R by the above assumptions. Then 
P satisfies the hypotheses (H) with fi replaced by /i(0 — fi{V^O/^V^- 

Remcirk 1. One example for such f is the following polynomial with only 
negative odd order terms 

fio = -aie — a^e^-' 

where ai > , i — 1, . . . ,p . 

Now we state the main result of this paper. 
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Theorem 2. Assume the Wiener process W satisfies assumption U^) . Then 
for any time T > and constant R > 0, there is positive constant C > , 
such that for any solution {u'' (t) , v"" (t)) of with initial value \ {uo, vq)\ < 

, there is a N -dimensional Wiener process W such that in distribution 

sup \u'{t) - y/euNiet) - ep7v(et)| < Ce^+ . (15) 

0<t<e-iT 

Here un solves ^ and solves the following stochastic differential equa- 
tions 

dtPN = A^pi, + VnUIMpn] + ^B{uN)dtW , (16) 
with pNif^) = and 

— 1 /"* 

foiuN) = lim - / foiuN + V*{s))ds 

t^oo t Jq 

and 

POO 

B{un) = 2/ ^ VNfo{uN + Tl*{t))-'PNTo{uN) 

® VNfoiuN + r]*{^)) -'PNfoiuN) dt, 
where (S> is the tensor product. Furthermore 

E sup \v'it)-^/^r|,{t)\<Ce^/^. (17) 

Now we introduce some scale transformations such that system ([Hll^dZl) 
is defined on time scale under the transformations. Introduce slow time 
t' = et and small fields = y/eu'''{x,t') and = y/ev'''{x,t') . Substituting 
and hereafter omitting primes, the coupled system (I6])-(I7!) is transformed 
to the following stochastic reaction diffusion equation with clearly separated 
time scales: 

dtu'' = d^^u' + VNriu' + v'), (18) 

dtv' = -d,.v' + QNr{u' + v') + ^QNdtW, (19) 
e ye 

with u'{0,t) = u'{Tr,t) = v'{0,t) = v'{TT,t) = for t > . The Wiener 
process is a rescaled version of ( fT2l) and with the same distribution. 

For convenience in the following we rewrite fll8p - (|T^ into the following 
form 



dtw' = A,w' + r{w') + J^AW , on < a; < tt . 



(20) 
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Here = [d^x , ^~^dxx) with zero Dirichlet boundary condition. 

In order to approximate solutions of ©^dZl), our basic idea is to pass 
limit e — s> in (fT8 l) -(fT9 l) to determine interesting structure in solutions, then 
study the deviation between the limit and solution. This first step is to 
give compact estimates, that is, tightness of solutions, as addressed in the 
following two sections. 

3 Stochastic convolution 

Start by considering the linear stochastic equation 

dtz' = A,z' + T.,dtW, z'{<d) = . 
Let Se(t) be the analytic semigroup generated by A^, then in a mild sense 

z'{t)= [ S,{t-s)J:jW{s). (21) 
Jo 

For any T > and 5 > , we give a uniform estimates for z{t), < t < T, 
in space Hq. We have 

Theorem 3. Assume ( 1731) . Then for any T > , q > , there is positive 
constant Cq{T) such that 

E sup \\z'{t)p < Cg{T). (22) 

o<i<r 

Proof. By using the stochastic factorization formula [7], for a G (0, 1/2) we 
have 

z\t) = ^-^^ / {t-sr-^S,{t-s)Y^{s)ds 



where 

ya{s)= I (s-r)-°5',(s-r)S,rfl^(r). (23) 
Jo 

Now for g > l/a , we have by the definition of A^ 



umi < c, {t-s)^ds- I \\s,{t-s)Y:{s)\\ids 

'-Jo -' -'0 

< c,{T) f\\Y^{s)\\lds 
Jo 

for some positive constant Cq{T). Then we have 

E sup \\z^{t)rs<c,{T)E r \\Y:{s)\\lds. 

0<t<T Jo 
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q/2 



q/2 



By the definition of A^, rewrite fl23|) as 

i=N+l ^0 

Then if g > 4, by the Bukholder-Davies-Gundy inequahty we have 
E||y;^(s)||^ = eII / (s-r)-"5,(s-r)S,diy(r) ' 

OO „g 

i=N+l •'^ 

OO 

i=l 

Therefore there is positive constant Cq(T), independent of e, such that 
E sup ||z^(t)||^<c,(T)E r \\Y:{s)rsds<C,{T). 

0<t<T Jo 

By the assumption ( |T3|) . taking 6 < 1 and the Young inequahty yields the 
result (122D for all g > . □ 

4 Tightness of solutions 

This section gives a tightness result by some a priori estimates for solutions 
to ffTSl) - ffTI?l) and the estimate of Theorem EJ 
Define w'^ = w"^ — z'^, then by 



dtw' = A,w' + f{w') (24) 

which is equivalent to 

with zero Dirichlet boundary condition on (0, tt) and w'^ = 11"^ + v'^ . 
Then we have for some positive constants C4 and C5 

—iwf = -\\ur-^~'\\vr + {nw^),w') 

^ II ~e||2 I e|2p , , ( \ e|2p-l , l/(2p)'\ \ \ 

< -\\w \\ - C2|u; |/2^(^) + CsTT + (^caltt; |/2p(/) + CsTT j |2;|^(2p)/(^) 
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Integrating with respect to time yields 

sup \w'{t)\^ + 2 [ \\w'{s)fds + 2ci [ \w'{s)\%,,j.ds 

0<t<T Jo Jo 

< \wo\' + 2c5 £ \z%s)\ll^y^j^ds + 2c,T 

< \wo\^ + 2c^f \\z\s)\\^Pds + 2c^T. (25) 

Jo 

On the other hand we have a positive constant C5 such that 

L(2p)'(/) 



Integrating with respect to time yields 

sup \\w\s)\\'' < ||M;oir + 2ci I sup \\w\T)\\^ds + 2 [ \w'{s)\%,,j. ds 

0<s<t Jo 0<T<s Jo 

+ 2c6^ \\z^{s)\\%^y^j^ds + 2ceT. 
Then by the Gronwall lemma and fl25|) . 

sup \\w'it)f<cjT(l + \\wof+ r \\z'is)f^ds] (26) 

0<t<T V ^0 / 

for some positive constant C7. Further for any integer m > 1 



dt 

Then 



sup ||u'''(s) 

0<s<t 



II 2m 

< Cm f sup ||«;^(r)f™rfs + C„ sup ||w)^(r)f— ^ f \w^{s)\%,,,.ds + 

Jo 0<T<s 0<T<T Jo 

rT 



a 



m 



Jo 



\\z^{s)f"''^ds + C^T+\\wof^. 
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By induction on m, we derive that 

sup \\w'it)f"'<CmT(l + \\wof'^+ r \\z'it)f'^''dt) (27) 

0<t<T ^ Jo ^ 

for some positive constant Cm- 

Now we show {C{w'')}e, the distribution of w*^, is tight in C(0, T; H). For 
this we need the following lemma by Simon [16J. 

Lemma 4. Assume E, Eq and Ei he Banach spaces such that Ei d Eq, 
the interpolation space (i?o, C E with 9 G (0, 1) and E d Eq with C 
and (E denoting continuous and compact embedding respectively. Suppose po, 
Pi E [1, oo] and T > 0, such that 

V is a bounded set in {0,T; Ei) 

and 

dV := {dv : V e V} is a bounded set in LP«(0, T; Eq). 
Here d denotes the distributional derivative. If 1 — 6 > 1/pe with 

I _i-e e 

Pe Po Pi ' 
then V is relatively compact in C(0,T; E). 

Now by the above lemma, noticing the estimate (|22l) and w"^ = w'^ — z"^, 
we draw the following result 

Theorem 5. Assume (C^j- For any T > , is tight in C{0, T; H) . 

5 Macroscopic reduction 

In this section we prove the main result. First Section 15.11 approximates 
the high modes by the Gaussian process ?7^. Then Section 15.21 derives an 
averaged approximation for the low modes, and the fluctuation is considered 
in Section 15.31 

5.1 Approximation for high modes 

We consider the high frequency dynamics of ffT^ . First for any fixed u G if at , 
v" satisfies 

dtv' = -d^.v' + Q^nu + v') + ^dtW. (28) 
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For any fo, f i G H]^, we have 

\v'it;vo)-v'it;v^)\'' < e-2--+i*/>o - ^^iT 

which means a unique stationary solution -y^ exists for any fixed u. In the 
following we determine the limit of in C(0,T;if^) as e ^ for any 
u E Hn . 

For this we scale time for equation (E]) by t' = tt which yields, upon 
omitting primes, 

dtV = -d,.V + -^dtW. (29) 

Here W is a rescaled version of the noise process in and with the same 
distribution. Then rjl{t) = rj*{t/e) is the unique stationary solution of ( !29l) . 
Moreover rjl is an exponential mixing Gaussian process with distribution 
;,=AA(0,a2(-A)-iQ2/2). 

Now we prove that for any u G Hjq , could be approximated by r/*" as 
e is small. Let = v"^ — t]^, then 

dtV' = -d^^^' + Q^fiu + t;^), ^^(0) = t;(0) - r/^(0). 
Multiplying on both sides of above equation in H yields 

||v..p<_f!^|v..|^,_i_i;.(„,„.)P. ,30) 

Then by the Gronwall lemma and fl7r|) . there is positive constant C such that 
for any t > 

g-a^+i(t-s)AE|/^(^ + v\s))\^ ds 



2Q!iV+l Jo 

< e'C{\\w{0)fP + E\r^Mn- 

Furthermore by 



for any T > , there is positive constant Ct such that 

E sup {t)-r]l{t)\<eCT {1^1^" + nvMD- 

0<t<T 

This proves (|T71) . 
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We end this subsection by giving an estimate on , m > , which 

is used in the fluctuation estimate in Section I5.3[ By (IHUjl we have 

dt € 2q;7v+1 

< -^^^\Vf- + eC^\f{u + v^)r 

for some positive constant Cm > . Then by the Holder inequahty, Gronwall 
lemma and (p7|) we have 

E\V'{t)\^"' <Cm, t>0. (31) 
5.2 Averaged equation 

In order to pass limit e ^ , we restrict our system into a small probability 
space. By the estimates in Section HI for any k > there is a compact set 
5« C C(0 , T ; H) such that 

¥{u' e B^} > 1 - K. 
Furthermore there is positive constant C^, such that 

\\u^{t)f<C^, te[0,T], 

for any G -Br. Now we introduce the probability space (fi^ , JF^ , P^) defined 
by 

= {uj en-.u' e B^}, j^^ = {snn^: s e J^} 

Then F{n <k. 

Now we restrict u & and introduce an auxiliary process. For any 
T > , partition the interval [0, T] into subintervals of length 6 = \fk . Then 
we construct processes such that for t G [^5, ik + 1)5), 

JkS 

u'iO)=uo, (32) 
dti'{t) = -d,.v'{t) + QNnu'{k6),v'{t)) + ^QNdtW{t), 

v'(kS) = v'{k5) . (33) 
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Then by the Ito formula for t G [k6, {k + 1)6), 
-^\v'(t)-d'(t)\^ 

< -^\v^(t) - v^it)\' + {nu\t) , v\t)) - nu\t) , nt))y{t) - m) 
+ {r{u\t) , v\t)) - r{u\k5) , v\t)) , v\t) - v\t)) 

< -^\v\t) - + e {c2\\w\t)\\^^ + C3) \u\t) - u\k5)\^ . 

By the choice of ^2^, there is Ct > , such that 

\u'{t) -u'{kS)\^ <CtS\ for te[k6,{k + 1)6). (34) 
Then by the Gronwall lemma, 

\v'{t) ~v'{t)\^ <Ct6\ te[0,T]. (35) 
In a mild sense for t G [k6, {k + 1)5) 

u'(t) = e^'^'-''^\'{kS) + f e^^'~'^VNr{u\s),v\s)) ds . 

JkS 

Then by the cubic property of / and smoothing property of e^*, noticing the 
choice of ^2^) we have for t G [k6, {k + 1)6) 

\u'{t) -u'{t)\ < C [ \v'{s)-v'{s)\ds + C' [ \u'{k6) -u'{s)\ds 

JkS JkS 

for some positive constant C . So by (IHSjl we have 

\u'{t) -u'{t)\ <Ct6 , tG[0,r]. (36) 
On the other hand, in a mild sense the solution of © is 

UM{t) = e^'uo + f e^^'~'^VNMuN{s)) ds . 
Jo 

Then, using [2;] to denote the largest integer less than or equal to z, 

{t)-UN{t)\ < f e^'^'-'^\VNr{u'{[s/6\6),ns))-VNT{u'{[s/6\6))\ds 
Jo 

+ [ e''^'-'^\VNnu^i[s/S\6))-V^Toiu'is))\ds 
Jo 

+ [ e^'^'-'^\rNTo{u'{s)) -VNTo{uN{s))\ds . 
Jo 
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Notice r] is independent of e, by the assumption H, and V^f" is continuous 
in e. Moreover the exponential mixing stationary measure n is independent 
of u, by the ergodic theorem 



/o(m7v) = hm - / fo{uN + 'n*{s))ds 



Then for any Ui , U2 E 



Hi, 



fo{uN + v)n{dv) . 



Hi 



[foiui,v) - fo{u2,v)]fi{dv] 



< 2cr 



< 2cn 



(■ui - U2){ul + ul + v^)fi{dv) 



Hi 



* l|2 



Ml - U2\ 



(37) 



which yields the continuity of VnIq '■ Hj^ — > H]\r . Then we have for t G [0, T] 



\u'{t) - UN{t)\ < Cj 



S+ \u^{s) — UN{s)\ds 



{31 



As 



\u'{t) - UNit)\ < \u'{t) - U{t)\ + \u{t) - UNit)\ , 

by the Gronwall lemma and ( l34l) . ( !36|) and ( !38|) we have for t G [0,T], 

|M^(t) -M7v(t)| < Ctv^. (39) 

Now by the arbitrariness of k, we complete the proof of the averaging 
approximation. And since T]^, is Gaussian with zero mean, we give 



VNfo{uN) = -CqVn{Uj^ + 3MArE?7 ) . 



(40) 



5.3 Fluctuation 

This subsection details the approximation of for small e. We study the 
deviation between u'' and u, which proves to be a Gaussian process. This 
shows that there are fluctuations in the slow modes. We follow an approach 
used previously [T71 |10l [18] . For this define the scaled difference 
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Then p*^ solves 

dtp' = d,,p' + ^ [VNriu', v') - Vn7{un)] , p'{0) = . 

However, here we just consider p*^, the solution of 

dtp' = d,,p' + [VnMu', v') - VnToM] , p'{0) = 0. 

By assumption (H) and estimate ( 1271) , E|p''(t) — p^it)] ^ as e ^ for any 
t > 0. 

Noticing the estimate flHTj) , by the Gronwall lemma for any T > , 

E sup |p^(t)r + E r \\p^it)fdt<CTil+\\wof). (41) 
o<t<r Jo 

In the mild sense we write 

p^{t) = 4 Te^^*-^) [VNfo{u'{r),v'{r)) - VNToiMr))] dr. 



Jo 

Then for any < s < t , by the property of e"^*, we have for some positive 
1 > 5 > 

\p'{t) - p'{s)\ < ^ te^^'-'^VNMu'{r),v'{r))~VNMuN{r))dr 



- r e^^'-'-^VNfoiu'ir),v^ir)) -VNfoMr))dr 
Jo 



< CT\t~sf^\VNfo{u%v')-VNfo{uN) 



L2(0,T;//jv) ■ 



By ( 14T1) and the estimates in Section H] 

E^JVNfo{u',v') -VnTo{un)\^,^,^t.h^^ < Cr(l + Ihof ). 

Then 

np'it)\c^O,T;H^)<CTil+\\Wof). (42) 

Here C\0,T] H^) is the H51der space with exponent 6. On the other hand, 
also by the property of e^*, we have for some positive constant Ct^o and for 
some 1 > a > 

\p'{t)\H'^ < f\t-s)-''/'\VNfo{u'{s),v'{s))~VNTo{uN{s))\ds 



'e 





< CT,a-^\PNfo{u% V') - ^Jv/o(«7v) |i2(o,T;Hjv) 



Macroscopic reduction for SRDEs 



16 



Then 

E sup \p'{t)\H^ < CtA^ + \\wo\\') . (43) 

0<t<T 

And by the compact embedding of C^(0, T; i/^)nC(0, T; H^) C C(0, T; Hn), 
the distribution of {p^}e is tight in C(0,T; ff). 
Spht = pI + p\ where each component satisfies, respectively, 

Qtp\ = QxxP\ + 4= [^iv/o(«iv, r;') - Pjv/o(Miv)] , Pi(0) = > 

= d..p\ + ^ [7^iv/o(^i^ ^1 - VNh{uN. rt)\ , ^2(0) = . 

Denote by v\ the probabihty measure of p\ induced on space C(0, T; -f/^iv)- 
And for 7 > denoted by UC^iH^.^ the space of all functions from iJiv 
to M which are uniformly continuous on together with all Frechet deriva- 
tives to order 7. In the following, for any h e UC'^{Hn,'R.) and un £ , 
denote by {h'{uN),-) '■ — > M the linear map defined by the first order 
Frechet derivatives of h and /i"(-ujv)(-) : ® Hn K the linear map 
defined by the second order Frechet derivatives of h . Then we have 

Lemma 6. Any limiting measure of vl, denote by , solves the following 
martingale problem on C(0, T; Hn): P°{pi(0) = 0} = 1 , 

h{p,{t))-h{p^m-j\h'{p,{T)),Ap,{T))dT-^j\r [h"{p^{T)){B{ur,))]dT 

is a P'^ -martingale for any h e UC^{Hn,'R). Here 



00 



B{un) = 2/ E[(Pjv/oK,#))-PAr/oK)) 

Jo 

® (Piv/oK,r;(0)) - VnTo{un))] dt. 
Proof. For any < s < t < 00 and h e UC°°{H) we have 

h{pl{t)) - h{p{{s)) 

{h'{pliT)),ApliT))dT 

^j\h\p\{T)),VNfo{uM{r)X{r)) - VmTo{um{t))) dr. 



[ 



+ 
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Rewrite the second term as 



1= j\h'{p{{T)),VNMuN{r),v'{r)) - VMTo{uN{r))) dr 

= 1= j\h'{pl{t)),VNMuN{r),v'{T)) -VNMuNir))) dr 

- i=^2*^''(pl(5))(p^/oK(r),r;^(T)) -P^^K(t))) 
<S)Apl{S)^ dSdr 

-\ jy^h'\p\{5)){vNh{uN{T),r)\T)) - VnTo{un{t)) 



®VnMun{S),v'{6)) -VNMuNm) dSdr 

Let {ci}^! be one eigenbasis of H, then 

h"ipl{S))(VNfo{uN{r),v'{r)) -VNTo{uN{r)) ^ 

N 

= d.Mp\{5)){{VNh{uN{T),v\T)) - VnToMS))} 

® {VNfoiuNiS),rj^iS)) - VnMun)}, O e^-^ . 

Here dij = de^dej where de^ is the directional derivative in direction e^. 
Denote by 

^W,r) = {{VNfo{uN{T),v'ir))-VNToMrm 

® {VNMuN{6),r]'{6)) - VnMM^))}: ei ® e,>. 

Then we have 

^3 = ff 9ijh{pl{S)){A'{S,r)ei,ej)dSdT 



X r) dXdSdr 
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e ,j JsJr 
^ ■ . J sJ s 



— L31 + L32 + L33 

where A^j{S,T) = A^j{6,T) — K[A^j{S,T)]. For our purpose, for any bounded 
continuous function $ on C{0,S]H), let $(-,ci;) = $(p|(-,u;)). Then by the 
exponential mixing of rj"^, 

|E[(L3i + L32)<I>]| -^Oase^O. 

Now we determine the limit of KAlj{6, t) d6 a.s e —>■ . For this introduce 

{Vj,foiu{r),r^^iT))-Vj,Toiuir))} 



Then 



E[A^^.(<5,r)-4.(5,r)] d5 



< 



E 



X (Wo(^(5),r/15)) - P7v/o(^(r),r/^(5)),e,)' 
By the assumption H we have 



d6 



\{VNMu{6),v'{S))-VNMu{T),v^{6)),e,)\ 
< 2[|m5)|p + h(r)f]K5)-t.(r)||e, 



and by dSZD 



|(Piv/o(w(r))-P^/o(w(5)),e,)| 



< 2 



l^'^)!!' + lk(T)f + E||r/f K<5)-^r)||e,|. 



Then also by the exponential mixing property of t]'^ 



l\[A^^{6,T)-A^^{6,T)]d6 



0, e ^ 



(44) 
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Now we put 



b^^{6-T) = E[{{VMfo{uN,Vm-VMfo{uN)} 



Then 



E[Al^{6,r)]=b^^ 



S-T 



Further, by the exponential mixing property, for any fixed 6 > t 



(5-r)/e 



Then, if e„ — as n ^ oo , z/''" P° , 

Jim E[L3$] = i rE^"(tr [h"{pr{r))B{uN)]<^)dT 



with B{un) = Bij{uN)ei ® cj . Similarly by the exponential mixing of r^" 

E[Li<l> + L2$] ^0 ase^O. 

By the tightness of in C{0,T; H), the sequence pi" has a limit process, 
denote by pi, in the weak sense. Then 



lim E 

n— >oo 



(/i'(p-(r)),M"M)<f^r 



E 



(/i'(pi(r)),M(T))$c^r 



and 



lim E[(/.(p-(t)) - /i(pr(.)))$] = E[(/.(pi(t)) - /i(pi(.)))$] . 

Jl^OO 

At last we have 

E^"[(%i)(i)-MPi(^)))*] 



E 



pO 



pO 



(/i'(pi(r)),MM)<^'^^r 

tr [h"{pi{r))BM]^dT 



(45) 



By an approximation argument we prove (H5l) holds for all /i G UC^{H). 
This completes the proof. □ 
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By fHUl) we have a more explicit expression of B{u]^) as 

/>oo 

B{un) = 2E / [P,v(3Miv(r?'-Er/2) + 3M^r7 + r7=^)] (46) 
Jo 

® [VN{^UN{rf{^) - ^rf) + 3M^r7(0) + 77=^(0))] dt . 

Then by the relation between weak solution to spdes and the martin- 
gale problem [TT], P° uniquely solves the martingale problem related to the 
following stochastic differential equation 

dtp^ = Ap^ + ^B{uN)dtW, pi(0) = 0, (47) 

where W{t) is A^-dimensional standard Wiener process, defined on a proba- 
bility space (OjjF, P) such that p\ converges weakly to pi in C{0,T; Hj^). 

By earlier results [17], P2 converges weakly to p2 in C{0,T; Hf^) and p2 
uniquely solves 

dtP2 = Ap2 + Vn[ToM{pi + P2)], P2(0) = 0. (48) 

Furthermore by foiu^) = —CoU%, we have fo{uj\f) = — 3coEf/^ . Then p^ 
converges weakly in C{0,T] Hjy) to p^ which uniquely solves the following 
A^-dimensional stochastic differential equation 

dtPN = ApN - ScoVn [uIpn + (Ef )p^] + ^/B{un) dtW , pn{0) = • 

6 Example of stochastic force in one mode 

This section applies the previous results to a simple case to see one example 
of how the noise forcing of high modes feeds into the dynamics of low modes. 
Further, we compare the result with that of the stochastic slow manifold 
model. For simplicity we assume the stochastic force acts just on the second 
spatial mode sin2x. 

6.1 Averaging and deviation 

Consider the following stochastic forced heat equation on the domain [0 , vr] 

dtW = d^^w + (1 + e-f)w -w^ + ay/edtW (49) 

with w{0, t) = w(7r, t) = . 7 is a real bifurcation parameter. The spatiotem- 
poral noise W is defined by ( IT2il with A2 = 1 and Aj = for i 7^ 2 , that is, 
W{x, t) = f32{t) sin2x . Then only the second spatial mode is forced by white 
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Figure 1: one realisation of the space-time dependence of the stochastic 
field w{x,t) for parameters 67 = y/ea = 1 . The stochastically forced sin2x 
mode interacts nonlinearly with the finite amplitude fundamental mode sinx . 
The numerics are based on finite differences with 15 points in space. 



noise. Figure [T] plots one realisation illustrating the nonlinear dynamics in- 
duced by the noise of strength ^/ea. In the spde fHUj) we incorporate the 
growth linear in w to counteract the dissipation on a finite domain so that 
we can control the clarity of the separation between fast and slow modes. 
We take Ai = dxx + 1 with Dirichlet boundary condition on [0,7r], then the 
eigenmodes ei{x) = sin(ix) corresponding to decay rates ai = P — 1 : giving 
the slow mode sin x ; and the fast modes sin ix for i > 2 . As the parameter €7 
crosses zero with no noise, a = 0, there is a deterministic bifurcation to a 
finite amplitude of the fundamental mode sin x . 

We consider the stochastic system (H9l) on long timescales of order . 
First notice that here f{w) = ew — , but all the analysis in Section [5] 
holds. Then decompose the field w'^{t) = y/eu''{t') + \/ev^{t') in the slow time 
t' = et . By Theorem [2] the fast mode v"^ is approximated by a stationary 
process rj which solves the following linear equation for small e > 

dt'T] = -dxxV + -^QnOvW . 

Decomposing rj = "^^riiei, then r^j = for z 7^ 2 and the scalar stationary 
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bifurcation parameter €7 

Figure 2: average equilibrium amplitude (squared) of the fundamental 
mode, a^, versus bifurcation parameter €7 for fixed noise amplitude ^/ea = 1 . 
The straight line fit is ^ +1.41e7 - 0.32 . 



process 772 then satisfies the following stochastic ordinary differential equation 

3 [T 
dm = — m dt' + crx - djS'^ . 
e V e 

The distribution of rj2 is the one dimensional normal distribution A/'(0 , |cr^) . 
Theorem [2] also asserts that the averaged equation for is 

dfU = •yu — Viu^ — ^^Viu sin^ 2x . (50) 

6 

Suppose u = A{t') sinx, by (I3UI) the amplitude A satisfies the Landau equa- 
tion 

— = 7 ]A A\ (51) 

dt' \' 4 J 4 ^ ^ 

Figure [2] plots the average amphtudes, a, of the fundamental mode obtained 
from numerical simulations. The figure shows that the amplitude of the fun- 
damental is depressed by the noise in the second component causing a delay 
in the bifurcation in the presence of noise. More extensive numerical simula- 
tions suggest that the mean amplitude depends upon noise and bifurcation 
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time t 

Figure 3: realisations of the amplitude of the fundamental mode a{t) versus 
time for various bifurcation parameters £7 and for noise amplitude ^/ea = 1 . 
At time zero the fundamental modes are started at the deterministic stable 
equilibrium for the corresponding €7. 



parameter (for < 67, ecr^ < 1) according to 

^ e(l.32 7 - 0.34 a^) + 6^(0.08 7^ + 0.03 7^2 + 0.02 a^). (52) 

The analytic Landau equation predicts equilibrium amplitude = ^•y — 
|cr^ which (after scaling by e) agrees remarkably well with these numerical 
estimates. 

However, there is also a significant stochastic component in the funda- 
mental mode, sinx, as shown by Figure [3] which plots the amplitude of the 
fundamental as a function of time for various bifurcation parameters €7. The 
stochastic fluctuations come from nonlinear interactions of the noise in the 
sin 2x mode. Now we calculate this deviation. Noting that f] = r]2 sin 2x , 
Lemma El asserts 

poo 

B{A) = 2E / [Vi{u + ri{s)f - Vi{u + ri)^] [Vi{u + ri{0)f -Vi{u + r]f] ds 
Jo 

POO 
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Then writing lim^^Q^u" — u)/y/e := pisinx, by 3Pi[(E?7^)pi] = cr^pi/4, the 
deviation pi solves the Ornstein-Uhlenbeck-like SDE 



9 . a' 



dpi= ~ - - -A' j p^dt + ^Ad(3 , pi(0) = 0, (53) 

where /3 is a standard real valued Brownian motion. For example, after the 
mean amplitude A reaches equilibrium, A = a/4(7 — a'^/A)/3, the SDE 
predicts fluctuations in pi with a standard deviation 

2 

(Ti ^ ^ = 0.1179 ea^ (54) 

Such fluctuations are seen in the numerical simulations of Figure [31 More 
extensive numerical simulations estimates the standard deviation of the fluc- 
tuations; Figure m plots this standard deviation against the noise amplitude^ 
A straight line fit to this data gives the standard deviation of the numeri- 
cally observed fluctuations as aa ~ 0.08 cr^. The theoretical prediction fl3^ 
scales the same with applied noise a, although the coefficient is about 30% 
different, and is similarly independent of the bifurcation parameter €7. Aver- 
aging and deviation together reasonably predict the dynamics of this example 
SPDE (|49l). 



6.2 Compare with the stochastic slow manifold 

Earlier work constructing stochastic slow manifolds of dissipative spdes ^1] 
is easily adapted to the example SPDE (HUj) . Recall we choose W = P2(t) sin 2x 
In terms of the amplitude a{t) of the fundamental mode sin a; , computer al- 
gebra readily derives that the stochastic slow manifold of the SPDE (H9l) is 

w = a sinx -|- — sin 3x + i/eo" sin 2x e~'^* * /32 

+ e^/%a sin 2x e"^* * e"^* i< $2 + ■ ■ ■ ■ (55) 

The history convolutions of the noise, e~^* * /?2 = j^_^^^^^^~^^dl32{s) that 
appear in the shape of this stochastic slow manifold empower us to eliminate 
such history integrals in the evolution except in the nonlinear interactions 
between noises; here simply 

d = e^a - - (/?2e~^* * /?2) H • 

^For larger stochastic forcing, ecr^ > 0.5 at this bifurcation parameter £7, the standard 
deviation Ua appears to plateau. 
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0.040 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 

noise ea^ 

Figure 4: standard deviation of the amplitude a{t) of the fundamental 
mode sin x versus noise ea^ for fixed bifurcation parameter e'j = 1 . These 
are estimated from long time versions of the simulations shown in Figure [3l 
The straight line fit predicts Ua ~ 0.07ecr^. 
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Analogously to the averaging and devaition theorems, analysis of Fokker- 
Planck equations [HI [Hj then asserts that the cannonical quadratic noise 
interaction term in this equation should be replaced by the sum of a mean 
drift and an effectively new independent noise process. Thus the evolution 
on this stochastic slow manifold is 



da 



a--a 



dt H —ea'^a d6 + ■ ■ ■ , (56) 

2^6 



where /3 is a real valued standard Brownian motion. The stochastic model ( l56l) 
is exactly the averaged equation ( 15T1) plus the deviation ( l53i) with a{t) = 
y/eA{et) +epi(et) . And the stochastic model (1521) predicts a stochastic equi- 
librium amplitude squared of about = |e7 — |eo"^ in agreement with the 
empirical fit (1521) to the numerical data of Figure [2] and other simulations. 

Now investigate the fluctuations about the stochastic equilibrium as seen 
in Figure [3] and measured in Figure HI As for the deviation equation ( |53ll . 
the stochastic slow model ( |56l) is approximately an Orstein-Uhlenbeck pro- 
cess in the vicinity of the finite amplitude stochastic equilibrium. Without 
elaborating the details, the form of fl3U|) then predicts fluctuations about the 
equilibrium have a standard deviation of ea'^/{Q\/2) in agreement with (IMl) 
and in moderate agreement with the fit of Figure H] to the numerical sim- 
ulations. The stochastic slow manifold model also predicts the stochastic 
dynamics of this example spde (1491) . The difference is that the stochastic 
slow manifold model is encapsulated in the one SDE (13U|) instead of being 
split into separate equations for the average and deviation. 
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